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Abstract 

We generalize the dynamical-mean field (DMFT) approximation by including into the DMFT 
equations some length scale £ via a momentum dependent "external" self-energy S^. This external 
self-energy describes non-local dynamical correlations induced by short-ranged collective SDW- 
like antiferromagnetic spin (or CDW-like charge) fluctuations. At high enough temperatures these 
fluctuations can be viewed as a quenched Gaussian random field with finite correlation length. 
This generalized DMFT+Sk approach is used for the numerical solution of the weakly doped one- 
band Hubbard model with repulsive Coulomb interaction on a square lattice with nearest and next 
nearest neighbour hopping. The effective single impurity problem in this generalized DMFT+Sk 
is solved by numerical renormalization group (NRG). Both types of strongly correlated metals, 
namely (i) doped Mott insulator and (ii) the case of bandwidth W < U (U — value of local 
Coulomb interaction) are considered. Densities of states, spectral functions and ARPES spectra 
calculated within DMFT+Sk show a pseudogap formation near the Fermi level of the quasiparticle 
band. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.+h, 74.72.-h 
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I. INTRODUCTION 



Among the numerous anomalies of the normal phase of high-temperature superconduc- 
tors the observation of a pseudogap in the electronic spectrum of underdoped copper oxides^ 
is especially interesting. Despite continuing discussions on the nature of the pseudogap, the 
preferable "scenario" for its formation is most likely based on the model of strong scattering 
of the charge carriers by short-ranged antiferromagnetic (AFM, SDW) spin fluctuations^. 
In momentum representation this scattering transfers momenta of the order of Q = (-, -) 
(a — lattice constant of two dimensional lattice). This leads to the formation of struc- 
tures in the one-particle spectrum, which are precursors of the changes in the spectra due 
to long-range AFM order (period doubling). As a result we obtain non-Fermi liquid like 
behaviour (dielectrization) of the spectral density in the vicinity of the so called "hot-spots" 
on the Fermi surface, appearing at intersections of the Fermi surface with antiferromagnetic 
Brillouin zone boundary (umklapp surface)- 

Within this spin-fluctuation scenario a simplified model of the pseudogap state was 
studied^** 1 under the assumption that the scattering by dynamic spin fluctuations can be 
reduced for high enough temperatures to a static Gaussian random field (quenched disor- 
der) of pseudogap fluctuations. These fluctuations are defined by a characteristic scattering 
vector from the vicinity of Q, with a width determined by the inverse correlation length 
of short-range order k = and by appropriate energy scale A (typically of the order of 
crossover temperature T* to the pseudogap stated). 

Undoped cuprates are antiferromagnetic Mott insulators with U ^> W (U — value of 
local Coulomb interaction, W — bandwidth of non-interacting band), so that correlation 
effects are very important. It is thus clear that the electronic properties of underdoped 
(and probably also optimally doped) cuprates are governed by strong electronic correlations 
too, so that these systems are typical strongly correlated metals. Two types of correlated 
metals can be distinguished: (i) the doped Mott insulator and (ii) the bandwidth controlled 
correlated metal W ~ U. Both types will be considered in this paper. 

A state of the art tool to describe such correlated systems is the dynamical mean-field the- 
ory (DMFT) 6 i 7 i 8 i 9 i 10 . The characteristic features of correlated systems within the DMFT are 
the formation of incoherent structures, the so-called Hubbard bands, split by the Coulomb 
interaction U, and a quasiparticle (conduction) band near the Fermi level dynamically gen- 



2 



erated by the local correlation a 6 i 7 i 8 i 9 i 10 . 

Unfortunately, the DMFT is not useful to the study the "antiferromagnetic" scenario of 
pseudogap formation in strongly correlated metals. This is due to the basic approximation 
of the DMFT, which amounts to the complete neglect of non-local dynamical correlation 
effects. 

Besides the extended DMFT—, which locally includes coupling to nonlocal dynamical 
fluctuations, a straightforward way to extend the DMFT are the so-called cluster mean-field 
theories^ 2 -. Two variants of this approach are the dynamical cluster approximation (DCA) 12 
and cellular DMFT (CDMFT) 1 -. In particular the DCA has been applied to study the low- 
energy properties of the Hubbard model, systematically including short- to medium ranged 
nonlocal correlations. Both improve on the cluster perturbation theory (CPT ) 14 ' 15 , a first 
attempts to use finite-size calculations to obtain approximate results for the thermodynamic 
limit. 

However, these approaches have certain drawbacks from both technical and interpreta- 
tion points of view. First, the effective quantum single impurity problem becomes rather 
complex. Thus, most computational methods available for the DMFT can be applied for the 
smallest clusters onl y 12 i 16 i 17 , i.e. include nearest-neighbor fluctuations only. For medium- to 
long-ranged correlations one is currently restricted to Quantum Monte-Carlo^ 8 -. Since for 
cluster problems again a sign problem arises, one is restricted to relatively small values of 
the local Coulomb interaction and high temperatures. Second, the interpretation of elec- 
tronic structures found has to be based on reliable input from other, typically approximate, 
complementary techniques. 

The aim of the present paper is to propose such a novel approach, which on the one hand 
retains the single- impurity description of the DMFT, viz a proper account for local corre- 
lations and the possibility to use very efficient impurity solvers like NR C 20 ' 21 ; on the other 
hand, we include non-local correlations on a non-perturbative model basis, which allows to 
control characteristic scales and also types of non-local fluctuations. This latter point allows 
for a systematical study of the influence of non-local fluctuations on the electronic proper- 
ties and in particular provide valuable hints on physical origin and possible interpretation 
of results found in e.g. more refined theoretical approaches. 

The paper is organized as follows: In section |H] we present a derivation of the self- 
consistent generalization we call DMFT+Sk which includes short-ranged dynamical corre- 



3 



lations to some extent. Section UTTI describes the construction of the k-dependent self-energy, 
and some computational details are presented in section IIV Al Results and a discussion are 
given in the sections II VI Then the paper is ended with summary section together with 
overview of related recent approaches and results on pseudogap issue. 



II. INTRODUCING LENGTH SCALE INTO DMFT: DMFT+S k APPROACH 

The basic shortcoming of the traditional DMFT approach 6 ' 7 ' 8 ^ 10 is the neglect of momen- 
tum dependence of the electron self-energy. This approximation in principle allows for an 
exact solution of correlated electron systems fully preserving the local part of the dynamics 
introduced by electronic correlations. To include non-local effects, while remaining within 
the usual "single impurity analogy", we propose the following procedure. To be definite, 
let us consider a standard one-band Hubbard model from now on. The extension to multi- 
orbital or multi-band models is straightforward. The major assumption of our approach 
is that the lattice and Matsubara "time" Fourier transformed of the single-particle Green 
function can be written as: 

GkM = -— * , v p r , u = nT(2n + l), (1) 

where T,(iu) is the local contribution to self-energy, surviving in the DMFT, while Ek(^) 
is some momentum dependent part. We suppose that this last contribution is due to either 
electron interactions with some "additional" collective modes or order parameter fluctua- 
tions, or may be due to similar non-local contributions within the Hubbard model itself. 

To avoid possible confusion we must stress that Hk(iu)) can in principle also contain local 
(momentum independent) contributions, which obviously vanish in the limit of infinite di- 
mensionality al — > oo and are not taken into account within DMFT. Due to this fact there is 
no double counting of diagrams within our approach to the Hubbard model. This question 
does not arise at all if we consider Ek(iu;) appearing due to some "additional" interaction. 
More important is that the assumed additive form of the self-energy £(zu;) + Ek(iu;) im- 
plicitly corresponds to neglect of possible interference of these local (DMFT) and non-local 
contributions. Furthermore, both contributions to the total self-energy Yi(iu) + Ek(iu;) in- 
dividually obeye causality by construction. Thus, the sum and finally the propagator (JTJ) 
constructed from it are causal, too. 
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The self-consistency equations of our generalized DMFT+Sk approach are formulated 
as follows: 

1. Start with some initial guess of local self-energy H(iu)), e.g. H(iu) = 0. 

2. Construct Ek(za;) within some (approximate) scheme, taking into account interactions 
with collective modes or order parameter fluctuations which in general can depend on 
£(ia>) and \x. 

3. Calculate the local Green function 

G«{iu) = ^Y, iuj + fl _ £(k) _ E(iw) _ Ek( iuj) • (2) 

4. Define the "Weiss field" 

^'M^H + G^M. (3) 

5. Using some "impurity solver" to calculate the single-particle Green function for the 
effective single Anderson impurity problem, defined by Grassmanian integral 

G d (r - T i) = ^-[ DcfDc ia c ia (r)4 a {T') exp(-S cff ) (4) 



n T ) n n( T ) , ( 5 ) 



with effective action for a fixed site ("single impurity") i 

5'cfr=- / dri j dT 2 c i a{ri)Qo 1 {T 1 -r 2 )ct(r 2 ) + / drUrii 
Jo Jo Jo 

Z c s = J Dcf a Dci a exp(—S e {[), and (3 = T~ l . This step produces a new set of values 

6. Define a new local self-energy 

sm^^h-g^m. (6) 

7. Using this self-energy as "initial" one in step 1, continue the procedure until (and if) 
convergency is reached to obtain 

Gaiiuj) = G d {iuj). (7) 

Eventually, we get the desired Green function in the form of (|T|). where H(iuj) and Sk(za;) 
are those appearing at the end of our iteration procedure. A more detailed derivation of 
this scheme within a diagrammatic approach is given in the Appendix EI 



III. CONSTRUCTION OF K DEPENDENT SELF ENERGY 



For the momentum dependent part of the single-particle self-energy we concentrate on 
the effects of scattering of electrons from collective short-range SDW-like antiferromagnetic 
spin (or CDW-like charge) fluctuations. To calculate Sj c (za;) for an electron moving in 
the quenched random field of (static) Gaussian spin (or charge) fluctuations with dominant 
scattering momentum transfers from the vicinity of some characteristic vector Q ("hot- 
spots" model 2 ), we use a slightly generalized version of the recursion procedure proposed 
in Refs . 4 i 5 i 22 which takes into account all Feynman diagrams describing the scattering of 
electrons by this random field. This becomes possible due to a remarkable property of our 
simplified version of "hot-spots" model that under certain conditions the contribution of an 
arbitrary diagram with intersecting interaction lines is actually equal to the contribution of 
some diagram of the same order without intersections of these line#^. Thus, in fact we can 
limit ourselves to consideration of only diagrams without intersecting interaction lines, taking 
the contribution of diagrams with intersections into account with the help of additional 
combinatorial factors, which are attributed to "initial" vertices or just interaction lines^ 2 -. 
As a result we obtain the following recursion relation (continuous fraction representation 2 ^) : 



iuj + /J, — H(iuj) — £ n (k) + inv n K — E„ +1 (ia;, k) 
Term E n (za;, k) of recurring sequence contains all contributions of diagrams with the number 

of interaction lines > n. Then 



is the sum of all diagrammatic contributions up to 2n-th order. Since the convergence of this 
recursion procedure for E n (iu;,k) is rather fast, one can take contribution for large enough 
n equal to zero and doing recursion backwards to n = 1 get desired physical self-energy^. 

The quantity A characterizes the energy scale and k = £ _1 is the inverse correlation length 
of short-range SDW (CDW) fluctuations, e n (k) = e(k + Q) and v n = \v£+q\ + \ v k+o\ f° r 
odd n while e n (k) = e(k) and v n = \v£\ + \ v^\ for even n. The velocity projections t>£ and 
are determined by usual momentum derivatives of the "bare" electronic energy dispersion 
e(k). Finally, s(n) represents a combinatorial factor with 




(9) 




(10) 
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for the case of commensurate charge (CDW-type) fluctuations with Q = (ir / 'a,n / 'a)&. For 
incommensurate CDW fluctuations 22 (when Q is not "locked" to the period of inverse lattice) 
one finds 

, s f 2±i for odd n 
s{n) = l 2 (11) 

| for even n. 

If we want to take into account the (Heisenberg) spin structure of interaction with spin 
fluctuations in "nearly antiferromagnetic Fermi-liquid" (spin-fermion (SF) model of Ref.—, 
SDW-type fluctuations), the combinatorics of diagrams becomes more complicated. Spin- 
conserving scattering processes obey commensurate combinatorics, while spin-flip scattering 
is described by diagrams of incommensurate type ("charged" random field in terms of Ref.— ). 
In this model the recursion relation for the single-particle Green function is again given by 
(JBJ, but the combinatorial factor s(n) now acquires the following forrn^: 

f for odd n 
s{n) = l 3 (12) 
I | for even n. 

Obviously, with this procedure we introduce an important length scale £ not present in 
standard DMFT. Physically this scale mimics the effect of short-range (SDW or CDW) 
correlations within fermionic "bath" surrounding the effective single Anderson impurity of 
the DMFT. We expect that such a length-scale will lead to a competition between local and 
non-local physics. 

An important aspect of the theory is that both parameters A and £ can in principle 
be calculated from the microscopic model at hand. For example, using the two-particle 
selfconsistent approach of Ref.^ 3 - with the approximations introduced in Refs4»^, one can 
derive within the standard Hubbard model the following microscopic expression for A: 

A = -U 2 ^— ^ < n iT > + < 71*1 > -2 < n^rin > = 



U = < n iT - riii >= 

= V^B^l < S? >, (13) 



where we consider only scattering from antiferromagnetic spin fluctuations. The different 
local quantities - spin fluctuation < Sf >, density n and double occupancy < n^riii > - 
can easily be calculated within the standard DMFT 9 . A detailed derivation of (|13|) and 
computational results for A obtained by DMFT using quantum Monte-Carlo (QMC) to 



solve the effective single impurity problem are presented in Appendix [Bj A corresponding 
microscopic expressions for the correlation length £ can also be derived within the two- 
particle self-consistent approach^. However, we expect those results for £ to be less reliable, 
because this approach is valid only for relatively small (or medium) values of U / 1. Thus, in 
the following we will consider both A and especially £ as some phenomenological parameters 
to be determined from experiments. 

IV. RESULTS AND DISCUSSION 
A. Computation details 

In the following, we want to discuss results for a standard one-band Hubbard model on 
a square lattice. With nearest (t) and next nearest it') neighbour hopping integrals the 
dispersion then reads 

e(k) = — 2t(cos k x a + cos k y a) — At' cos k x a cos k y a , (14) 

where a is the lattice constant. The correlations are introduced by a repulsive local two- 
particle interaction U. We choose as energy scale the nearest neighbour hopping integral t 
and as length scale the lattice constant a. 

For a square lattice the bare bandwidth is W = 8t. To study a strongly correlated metallic 
state obtained as doped Mott insulator we use U = 40i as value for the Coulomb interaction 
and a filling n = 0.8 (hole doping). The particular choice of the latter value for U is 
motivated by two aspects. First, this value of U leads to an insulating DMFT+Ek solution 
at half-filling. Second, estimations of U for stoichiometric La2Cu04 (high-Tc prototype 
compound) based on constrained LDA 2 » calculations typically give U of the order of 10 eV^i, 
which corresponds to 40t to our choice of parameters. The correlated metal in the case of 
W > U is realized via U — 4t - a value used in various theoretical papers discussing the 
pseudogap state - and two fillings: half-filling (n = 1.0) and n = 0.8 (hole doping). As 
typical values for A we choose A = t and A = 2t (actually as approximate limiting values 
- cf. Appendix |B} and for the correlation length £ = 2a and £ = 10a (motivated mainly 
by experimental data for cuprates^). 

The DMFT maps the lattice problem onto an effective, self-consistent single impurity 
defined by Eqs. 0-©. In our work we employ as "impurity solvers" two reliable nu- 
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merically exact methods — quantum Monte-Carlo (QMC) 18 and numerical renormalization 
group (NRG) 20 i 21 . Calculations were done for the case t' — and t'/t=-0A (more or less 
typical for cuprates) at two different temperatures T = 0.088t and T = 0.356t (for NRG 
computations) 40 . QMC computations of double occupancies as functions of filling were done 
at temperatures T = O.lt and T = OAt— . 

Below we present results only for most typical dependences and parameters, more data 
and figures can be found in Ref.— . 

B. Generalized DMFT+Ek approach: densities of states 

Let us start the discussion of our results obtained within our generalized DMFT+Sk 
approach with the densities of states (DOS) for the case of small (relative to bandwidth) 
Coulomb interaction U = At with and without pseudogap fluctuations. As already discussed 
in the introduction, the characteristic feature of the strongly correlated metallic state is the 
coexistence of lower and upper Hubbard bands split by the value of U with a quasiparticle 
peak at the Fermi level. Since at half-filling the bare DOS of the square lattice has a 
Van-Hove singularity at the Fermi level (f = 0) or close to it (in case of t'/t = —0.4) one 
cannot treat a peak on the Fermi level simply as a quasiparticle peak. In fact, there are two 
contributions to this peak from (i) the quasiparticle peak appearing in strongly correlated 
metals due to many-body effects and (ii) the smoothed Van-Hove singularity from the bare 
DOS 42 . In Figs, ^and El we show the corresponding DMFT(NRG) DOS without pseudogap 
fluctuations as black lines for both bare dispersions t'/t = —0.4 (left panels) and for t' — 
(right panels) for two different temperatures T = 0.356t (lower panels) and T = 0.088t 
(upper panels) with fillings n = 1.0 and n = 0.8 respectively. The remaining curves in 
Figs. Qand 121 represent results for the DOS with non-local fluctuations switched on with the 
fluctuation amplitude A = 2t. For all sets of parameters one can see that the introduction 
of non-local fluctuations into the calculation leads to the formation of pseudogap in the 
quasiparticle peak. 

The behaviour of the pseudogaps in the DOS has some common features. For example, 
for t'=0 at half-filling (Fig. right column) we find that the pseudogap is most pronounced. 
For n = 0.8 (Fig. [21 right column) the picture is almost the same but slightly asymmetric. 
The width of the pseudogap (the distance between peaks closest to Fermi level) appears to 
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be of the order of ~ 2A here. Decreasing the value of A from 2t to t leads to a pseugogap 
that is correspondingly twice smaller and in addition more shallow (see Ref.—). When one 
uses the combinatorial factors corresponding to the spin-fermion model fEq.(|12j)). we find 
that the pseudogap becomes more pronounced than in the case of commensurate charge 
fluctuations (combinatorial factors of Eq. (jlljl ). The influence of the correlation length £ 
can be seen is also as expected. Changing from = 0.1 to = 0.5, i.e. decreasing the 
range of the non-local fluctuations, slightly washes out the pseudogap. Also, increasing the 
temperature from T = 0.088t to T = 0.356/f: leads to a general broadening of the structures 
in the DOS. These observations remain at least qualitatively valid for t'/t = —0.4 (Figs, 
n and |2l left columns) with an additional asymmetry due to the next-nearest neighbour 
hopping. Noteworthy is however the fact that for t'/t = —0.4 and £ _1 = 0.5 the pseudogap 
has almost disappeared for the temperatures studied here. Also very remarkable point is 
the similarity of the results obtained with the generalized DMFT+Ek approach with U = At 
(smaller than the bandwidth W) to those obtained earlier without Hubbard-like Coulomb 
interactions^. 

Let us now consider the case of a doped Mott insulator. The model parameters are 
t'/t = —0.4 with filling n = 0.8, but the Coulomb interaction strength is now set to U — 
40i. Characteristic features of the DOS for such a strongly correlated metal are a strong 
separation of lower and upper Hubbard bands and a Fermi level crossing by the lower 
Hubbard band (for non-half-filled case). Without non-local fluctuations the quasi-particle 
peak is again formed at the Fermi level; but now the upper Hubbard band is far to the 
right and does not touch the quasiparticle peak (as it was for the case of small Coulomb 
interactions). DOS without non-local fluctuations are again presented as black lines in Fig. 01 
Results for the case if — are presented elsewhere^. 

With rather strong non-local fluctuations A = 2t, a pseudogap appears in the middle of 
quasiparticle peak. In addition we observe that the lower Hubbard band is slightly broadened 
by fluctuation effects. Qualitative behaviour of the pseudogap anomalies is again similar to 
those described above for the case of U = 4t, e.g. a decrease of £ makes the pseudogap less 
pronounced, reducing A from A = 2t to A = t narrows of the pseudogap and also makes 
it more shallow etc. (see Ref.—). Note that for the doped Mott-insulator we find that the 
pseudogap is remarkably more pronounced for the SDW-like fluctuations than for CDW-like 
fluctuations. 
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There are, however, obvious differences to the case with U = At. For example, the width 
of the pseudogap appears to be much smaller than 2 A, beeing of the order of A/2 instead 
(see Fig. OJ. This effect we attribute to the fact that the quasiparticle peak itself is actually 
strongly narrowed now by local correlations. 



C. Generalized DMFT+Sk approach: spectral functions A(uj,\s) 

In the previous subsections we discussed the densities of states obtained self-consistently 
by the DMFT+Sk approach. Once we get a self-consistent solution of the DMFT+Ek 
equations with non-local fluctuations we can of course also compute the spectral functions 
A(u,k) 

A(u, k) = --Im— — 1 (15) 

7r uj + n — e[k) — h(u) — h^u) 

where self-energy £(o>) and chemical potential /i are calculated self-consistently as described 
in Sec. |H] To plot A(u,\t) we choose k-points along the "bare" Fermi surfaces for different 
types of lattice spectra and filling n = 0.8. In Fig. 0]one can see corresponding shapes of 
these "bare" Fermi surfaces (presented are only 1/8-th of the Fermi surfaces within the first 
quadrant of the first Brillouin zone). 

A first natural quantity to inspect is the self-energy S(k, uj + i8), shown in Fig. El for 
t'/t = —0.4, n = 0.8 and U = At (left column) and U = AOt (right column). As repre- 
sentative k-points we chose the centre of the first Brillouin zone (r), the "hot-spot" and 
"cold-spot" (point "B" in Fig. UJ. The results were obtained with NRG at a temperature 
T = 0.0881 The structures for U = At are rather broad, but reveal after a closer inspection 
features similar to the case U = AOt. For the latter, the behaviour at T and "B" is very 
different from the structures at the "hot-spot" . Namely, while for the former two k-points 
ImS(k, uj + i5) shows a nice parabolic maximum at the Fermi energy, the latter develops a 
minimum instead. Such a structure in the self-energy will result in a rather evident (pseudo) 
gap in the spectral function at this k-point and weaker pseudogap behaviour in the DOS. 
Its appearance is obviously due to the presence of the spin-fluctuations at the "hot-spot". 
Note that similar features have been observed in numerically expensive cluster mean-field 
calculations 27 , too, with an interpretation as spin-fluctuation induced based on physical ex- 
pectations. Our calculations, obtained at a minimum numerical expense, indeed show, that 
including short-ranged fluctuations will precisely produce these non Fermi-liquid structures 
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in the one-particle self-energy. This behaviour is quite typical for the problem and was ob- 
served by other groups using different method o 16 ' 28 ' 29 ' 30 . In several works midgap peak in the 
pseudogap was obtained with explanation of its origin by particular shape of the self-energy 
close to the Fermi leve r 28 ' 29 ' 31 . 

In the following we concentrate mainly on the case U = At and filling n = 0.8 (Fermi 
surface of Fig. Ufa)). The corresponding spectral functions A(u,k) are depicted in Fig. El 
When if ft = —0.4 (upper row), the spectral function close to the Brillouin zone diagonal 
(point B) has the typical Fermi-liquid behaviour, consisting of a rather sharp peak close to 
the Fermi level. In the case of SDW-like fluctuations this peak is shifted down in energy 
by about — 0.5t (left upper corner). In the vicinity of the "hot-spot" the shape of A(u,\s) 
is completely modified. Now A(u, k) becomes double-peaked and non-Fermi-liquid-like. 
Directly at the "hot-spot", A(u,k) for SDW-like fluctuations has two equally intensive 
peaks situated symmetrically around the Fermi level and split from each other by ~ 1.5 A 
Refs.^. For commensurate CDW-like fluctuations the spectral function in the "hot-spot" 
region has one broad peak centred at the Fermi level with width ~ A. Such a merging 
of the two peaks at the "hot-spot" for commensurate fluctuations was previously observed 
in Ref.— . However close to point A this type of fluctuations also produces a double-peak 
structure in the spectral function. 

Spectral functions for the case of U = At at half-filling (n = 1) and for if ft = —0.4 are 
similar to those just discussed for n = 0.8. However, the pseugogap is more pronounced in 
this case and remains open everywhere close to the umklapp surface for SDW fluctuations^. 

In the lower panel of Fig. El we show spectral functions for 20% hole doping (n = 0.8) and 
the case of t' = (Fermi surface from Fig. Efb)). Since the Fermi surface now is close to 
the umklapp surface, the pseudogap anomalies are rather strong and almost non-dispersive 
along the Fermi surface. At half-filling for if = the Fermi surface actually coincides with 
umklapp surface (in case of perfect "nesting" whole Fermi surface is the "hot-region"). The 
spectral functions are now symmetric around the Fermi level. For SDW-like fluctuations 
there are two peaks split by ~ 1.5A. Again, CDW-like fluctuations give just one peak 
centred at the Fermi level with width ~ A. 

For the case of a doped Mott insulator (U = AOt, n = 0.8), the spectral functions obtained 
by the DMFT+E k approach are presented in Fig. [71 Qualitatively, the shapes of these 
spectral functions are similar to those shown in Fig. El As was pointed out above, the strong 
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Coulomb correlations lead to a narrowing of the quasiparticle peak and a corresponding 
decrease of the pseudogap width. As is evident from Fig. [7| the structures connected to the 
pseudogap are now spread in an energy interval ~ t, while for U = At they are restricted 
to an interval ~ At instead. One should also note that in contrast to U — At the spectral 
functions are now about four times less intensive, because part of the spectral weight is 
transferred to the upper Hubbard band located at about AOt and well separated from the 
quasiparticle peak now. 

Using another quite common choice of k-points we can compute A(cu, k) along high- 
symmetry directions in the first Brillouin zone: r(0, 0) — X(7r, 0) — M(n, it) — r(0, 0). The 
spectral functions for these k-points are collected in Fig. |H] for the case of SDW-like 
fluctuations. Characteristic curves for doped Mott insulator are presented in Ref . For all 
sets of parameters one can see a characteristic double - peak pseudogap structure close to 
the X point. In the middle of M — T direction (so called "nodal" point) one can see the 
reminiscence of AFM gap which has its biggest value here in case of perfect antiferromagnetic 
ordering. Also in the nodal point "kink" -like behaviour is observed caused by interactions 
between correlated electrons with short-range pseudogap fluctuations. A change of the 
filling leads mainly to a rigid shift of spectral functions with respect to the Fermi level. 

With the spectral functions we are now of course in a position to calculate angle resolved 
photoemission spectra (ARPES), which is the most direct experimental way to observe 
pseudogap in real compounds. For that purpose, we only need to multiply our results for 
the spectral functions with the Fermi function at temperature T = 0.088t. Typical example 
of the resulting DMFT+Sk ARPES spectra are presented in Fig. El More figures of ARPES- 
like results obtained within the DMFT+Sk approach for a variety of parameters can be found 
in Ref.— . One should note that for t'/t = —0.4 (upper panel of Fig. EJ) as k goes from point 
"A" to point "B" the peak situated slightly below the Fermi level changes its position and 
moves down in energy. Simultaneously it becomes more broad and less intensive. The dotted 
line guides the motion of the peak maximum. Also at the "hot-spot" and further to point 
"B" one can see some signs of the double-peak structure. Such behaviour of the peak in the 
ARPES is rather reminiscent of those observed experimentally in underdoped cuprate o 2 ! 4 ' 32 . 
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V. CONCLUSION 



To summarize, we propose a generalized DMFT+Sk approach, which is meant to take 
into account the important effects of non-local correlations (in principle of any type) in 
addition to the (essentially exact) treatment of local dynamical correaltions by the DMFT. 
In the standard DMFT the "bath" surrounding the effective single Anderson impurity is 
spatially uniform since the DMFT self-energy is only energy-dependent. The main idea 
of our extension is to introduce non-local correlations through the "bath", i.e. to make 
it spatially non-uniform, while keeping standard DMFT self-consistency equations. Such a 
generalization of the DMFT allows to supplement it with a k-dependent self-energy E(k, ui). 
It in turn opens the possibility to access the physics of low-dimensional strongly correlated 
systems, where different types of spatial fluctuations (e.g. of some order parameter) become 
important, in a non-perturbative way at least with respect to the important local dynamical 
correlations. However, we must stress that our procedure in no way introduces any kind of 
systematic l/(i-expansion, being only a qualitative method to include a length scale into 
DMFT. Nevertheless we believe that such a technique can give valuable insight into the 
physical processes leading to correlation induced /c-dependent structures in single-particle 
properties. 

In this work we model such effects for the two-dimensional Hubbard model by incorporat- 
ing into the "bath" scattering of fermions from non-local collective SDW-like antiferromag- 
netic spin (or CDW-like charge) short-range fluctuations. The corresponding k-dependent 
self-energy E(k, to) is obtained from a non-perturbative iterative scheme^. Such choice of 
the S(k, to) allows to address the problem of pseudogap formation in the strongly correlated 
metallic state. We showed evidence that the pseudogap appears at the Fermi level within 
the quasiparticle peak, introducing a new small energy scale of the order of psedogap poten- 
tial value A in the DOS and more pronounced in spectral functions A(u, k). Let us stress, 
that our generalization of DMFT leads to non-trivial and in our opinion physically sensible 
k-dependence of spectral functions. Is is significant that this particular choice of S(k, uj)r^ 
does not cause difficulties to "double counting" problem within our combined DMFT+Ek 
approach. Also, the combination of diagrammatically correct techniques like DMF T 6 ' 7 ' 8 ' 9 ' 10 
and the non-local self-energy ansatz of Refs4»£ preserves the correct analytical properties 
of the combined self-energy Y,(iuj) + T,y{iu), as well as of the corresponding one-electron 
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propagator (JTJ). 

Of course our pseudogap observations are not entirely new. Similar results about 
pseudogap formation in the 2d Hubbard model were already obtained within cluster 
DMFT extensions, i.e. the dynamical cluster approximation (DCA) 12 ' 27 and cellular DMFT 
(CDMFT) 16 i 17 , CP T 14 ' 15 ' 33 and two interacting Hubbard sites selfconsistently embedded in 
a bath 2 ^ However, these methods have generic restrictions concerning the size of the clus- 
ter, temperatures or filling accessible and, in case of the QMC, values of the local Coulomb 
energy. Recently, also the EDMFT was applied to demonstrate pseudogap formation in the 
DOS due to dynamic Coulomb correlations 34 . Note, however, that within the EDMFT there 
is no way to obtain a k-dependence in spectral functions beyond that originating from the 
bare electronic energy dispersion. Important progress was also made with weak coupling 
approaches for the Hubbard model 35 and functional renormalization group 29 ' 30 . In several 
papers pseudogap formation was described in the framework of the t-J model^. A more 
general scheme for the inclusion of non-local corrections was also formulated within the so 
called GW extension to the DMFT^ 7 -^. 

While at a first glance the introduction of additional phenomenological parameters (cor- 
relation length £, and A) through the definition of S(k, u) seems to be a step back with 
respect the methods outlined above, it actually opens the possibility to systematically distin- 
guish between different types of nonlocal fluctuations and their effects and helps to analyze 
experimental or theoretical data obtained within more advanced schemes in terms of intu- 
itive physical pictures. Note, however, that in principle even the paramters £ and A can be 
calculated from the original model^ 3 -, too. 

An essential advantage of the proposed combination of two non-perturbative methods 
(DMFT and S(k, u) from Refs4^) removes the restrictions on model parameters in e.g. 
cluster mean-field theories. Our scheme works for any Coulomb interaction strength U, 
pseudogap strength A, correlation length £, filling n and bare electron dispersion e(k) on 
a 2d square lattice for any set of k-points. Although we presented only high-temperature 
data in this paper, the possibility to use Wilson's NRG to solve the effective impurity model 
also opens the possibiltiy to study properties at T = 0, which is currently impossible within 
the DCA or CDMFT for larger clusters. Moreover, the DMFT+Sk approach can be easlily 
generalized to orbital degrees of freedom, phonons, impurities, etc. 

As a further application of our generalized DMFT+Sk we would like to bring readers 
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attention to Ref^, dealing with the problem of Fermi surface destruction in High-T c com- 
pounds because of pseudogap fluctuations. 
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APPENDIX A: DERIVATION OF GENERALIZED DMFT+S k APPROACH 



In this appendix we present a derivation of the generalized DMFT+Ek scheme for the 
Hubbard model 

H = - *ii4r c *r + U n i\ n ih ( Al ) 
ij,cr i 

using a diagrammatic approach. The single-particle Green function in Matsubara represen- 
tation is as usual given by 

G ^ = -T ( L w (A2) 

ilu + fi — e(k) — h(iu, k) 

To establish the standard DMFT one invokes the limit of infinite dimensions d — > oo. In 
this limit only local contributions to electron self-energy survive^, i.e. E^- — > 5yEjj or, in 
reciprocal space, E(iu;,k) — > H(iu). 

In Fig. El we show examples of "skeleton" diagrams for the local self - energy, contribut- 
ing in the limit of d — > oo. The complete series of these and similar diagrams defines the 
local self - energy as a functional of the local Green function 

E = F[Gn] , (A3) 

where 

G*M = — * w . r - (A4) 

N ilj + fi — e{K) — H[zu) 

One then defines the "Weiss field" 

^M^M + G^M (A5) 

which is used to set up the effective single impurity problem with an effective action given by 
(JHJ). Via Dyson's equation the Green function (jlj) for this effective single impurity problem 
can be written as 

1 

and the "skeleton" diagrams for self-energy E^ are just the same as shown in Fig. with 
the replacement Gu — * Gd- Thus we get 



G dM = ^7—. (A6) 



S d = F[G d ], (A7) 



where F is the same functional as in (|A3|) . The two equations ()A6|) and (|A7|) define both 
Gd and E d for a given "Weiss field" Q . On the other hand, for the local E and Gu of the 
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initial (Hubbard) problem we have precisely the same pair of equations, viz (jA3|) and (|A5|) . 
and Go i n both problems is just the same, so that 

£ = Gu = G d . (A8) 

Thus, the task of finding the local self-energy of the (d — > oo) Hubbard model is eventually 
reduced to the calculation of the self-energy of an effective quantum single impurity problem 
defined by effective action of Eq. (JHJ). 

Consider now non - local contribution to the self - energy. If we neglect interference 
between local and non-local contributions (as given e.g. by the diagram shown in Fig lTlT b)). 
the full self-energy is approximately determined by the sum of these two contributions. 
"Skeleton" diagrams for the non-local part of the self-energy, Ek(iu;), are then those shown 
in Fig. rrTT a). where the full line denotes the Green function of Eq. while broken 
lines denote the interaction with static Gaussian spin (charge) fluctuations. These diagrams 
are just absent within the standard DMFT (as any contribution from Ornstein - Zernike 
type fluctuations vanish for d — > oo), and no double counting problems arise at all. 

The local contribution to the self-energy is again defined by the functional (jA3|) via the 
local Green function Gu, which is now given by (0). Introducing again a "Weiss field" via 
(IA5jl and repeating all previous arguments, we again reduce the task of finding the local part 
of the self-energy to the solution of an "single impurity" problem with an effective action 

©• 

To determine the non-local contribution H^iu) we first introduce 

£okM = r-rrr- \ = -■ — , 7T\ — w~T ( A9 ) 

G k \iuj) + H^iuj) ioj + /i — e{k) — h(tuj) 

as the "bare" Green function for electron scattering by static Gaussian spin (charge) fluctua- 
tions. The assumed static nature of these fluctuations allows to use the method of Refs . 4 ^ 122 
and the calculation of the non-local part of the self-energy H^iu) reduces to the recursion 
procedure defined by Eqs. © and ©• The choice of the "bare" Green function Eq. (|A9|) 
guarantees that the Green function "dressed" by fluctuations G^^cu) = G ^(i^j) — £k(^), 
which enters into the "skeleton" diagrams for ^(iu), just coincides with the full Green 
functions G^iu). 

Thus we obtain a fully self-consistent scheme to calculate both local (due to strong single- 
site correlations) and non-local (due to short-range fluctuations) contributions to electron 
self-energy. 
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APPENDIX B: A IN THE HUBBARD MODEL. 



In this Appendix we derive the explicit microscopic expression for pseudogap amplitude 
A given in (|13|). Within the two-particle self-consistent approach of Ref.— , valid for medium 
values of U, and neglecting charge fluctuations, we can write down an expression for the 
electron self-energy of the form used in (0), with 

E a {iu) = Un_ a (Bl) 

as the lowest order local contribution due to the on-site Hubbard interaction, surviving in the 
limit of d —>■ oo, and exactly accounted for in DMFT (with all higher-order contributions). 
Non-local contribution to the self-energy (vanishing for d — > oo and not accounted within 
DMFT) due to interaction with spin-fluctuations then leads to the expression 

U T . 

S fc(H = T at Y Y U "pX v {fh v m )G (k + q, iuj + iv m ) , (B2) 



4 N 

m q 

where 

U sp = g n (0)U, g n (0) = < > (B3) 

with < r? a >=< n a > and < >=< >= \n in the paramagnetic phase. For the 
dynamic spin susceptibility Xsp(q> v m) we use the standard Ornstein-Zernike forrn^ 3 ., similar 
to that used in spin-fermion model^, which describes enhanced scattering with momenta 
transfer close to antiferromagnetic vector Q = (n/a, n/a). With these approximations, we 
can write down the following expression for the non-local contribution to the self-energy^ 5 -: 

1 



1 T . 



4 N ^— ' ^—f ^ v "' m> ico + ip m + f ,- e (k + q) 



\ UU n>Tf Y Y Xs p^ Vm ^ Y 



N ^— ' ^— ' ^— ' icu + a — e(k + q) 
A 2 V5(q): 



q 

A 2 



V - e(p + Q) + «Op+ql + l^p +Q l)^sign^' 
Here we have introduced the static form factor 5 



(B4) 



m (a, - Q,) 2 + ^ 2 (g„ - Q,) 2 + 1~ 2 (B5) 
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and the squared pseudogap amplitude 



A 2 = -^UU sp — Xs p^ Vm ) 



m q 



1 



UU sp [< > + < n i]r > -2 < n iT n^ >] 



4 



\uu t 



sp 



(B6) 



where we have used the exact sum-rule for the susceptibility 4,23 . Taking into account (|B3|) 
we immediately obtain (fT3|). 

Actually, the approximations made in (JB4j) and (JB5j) allow for an exact summation of the 
whole Feynman series for electron interaction with spin-fluctuations, replaced by the static 
Gaussian random field. Thus generalizing the one-loop approximation (|B4j) eventually leads 
to the basic recursion procedure given in 0, © Refs.- 4 ^. 

Using the DMFT(QMC) approach we computed occupancies < >, < > and double 
occupancies < n^n^ > required to calculate the pseudogap amplitude A of Eq. (jB6|) In Fig. 
El the corresponding values of A are presented. One can see that A grows when the filling 
goes to n = 1. While U approaches 8t (the value of the bandwidth for a square lattice) A as 
a function of n grows monotonically. When U becomes larger than W = 8t (when a metal- 
insulator transition occurs) one can see a local minimum for n = 0.9, which becomes more 
pronounced with further increase of U. For t'/t = —0.4 and both temperatures the scatter 
of A values is smaller than for the case of t' — 0. Also A has a rather weak temperature 
dependence. All values of A lie in the interval ~ 0.75t-h2t. Therefore, for our computations 
we took only two characteristic values of A = t and A = 2t. 
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FIG. 1: (Color online) Comparison of DOS obtained from DMFT(NRG)+£k calculations for 
different combinatorical factors (SF — spin-fermion model, commensurate), inverse correlation 
lengths (£ _1 ) in units of the lattice constant, temperatures (T) and the value of pseudogap potential 
A = 2t. Left column corresponds to t'/t = —0.4, right column to t' = 0. In all graphs the Coulomb 
interaction is U = At and n = 1. The Fermi level corresponds to zero. 
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FIG. 2: (Color online) Comparison of DOS obtained from DMFT(NRG)+Sk calculations for a 
filling n = 0.8, other parameters as in Fig. ^ 
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FIG. 3: (Color online) Comparison of DOS obtained from DMFT(NRG)+£k calculations 
t'/t = -0.4, T = 0.088t, U = 40i, A = 2t and filling n = 0.8. 




FIG. 4: 1/8-th of the bare Fermi surfaces for the occupancy n = 0.8 and different combinations 
(t, t') used for the calculation of spectral functions A(\c,u). The diagonal line corresponds to the 
umklapp surface. The full circle marks the so-called "hot-spot" . 
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FIG. 5: (Color online) Real (dashed line) and imaginary (full line) parts of the self-energy E(k, u) 
for t/t' = —0.4, U = At (left column) and U = 40t (right column) for characteristic k-points: T, 
"hot-spot" (see Fig. 0J) and "cold-spot" (point "B" in Fig. gj). For all graphs the filling is n = 0.8, 
temperature T = 0.088i, inverse correlation length £ _1 = 0.1, value of pseudogap potential A = 2t, 
and SF combinatorics. 
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FIG. 6: Spectral functions A(k,uj) obtained from the DMFT(NRG)+£k calculations along the 
directions shown in Fig. |1J Model parameters were chosen as U = At, n = 0.8, A = 2t, £ _1 = 0.1 
and temperature T = 0.088t. The "hot-spot" k-point is marked as fat dashed line. The Fermi 
level corresponds to zero. 
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FIG. 7: Spectral functions A(k, w) obtained from the DMFT(NRG)+$]k calculations for U 
other parameters as in Fig. HO 
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FIG. 8: Spectral functions A(k,u>) obtained from the DMFT(NRG)+S k calculations along high- 
symmetry directions of first Brillouin zone r(0, 0) — X(ir, 0) — M(ir, ir) — T(0, 0), SF combinatorics 
(left row) and commensurate combinatorics (right column). Other parameters are U = 4i, n = 0.8, 
A = 2t, ^ 1 = 0.1 and temperature T = 0.088t. The Fermi level corresponds to zero. 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




-3 -2 -10 -2-10 1 

Energy (units of t) 

FIG. 9: ARPES spectra simulated by multiplication of the spectral functions obtained from 
DMFT(NRG)+S k calculations for U = At and n = 0.8 Fig. El with Fermi function at T = 0.088* 
plotted along the lines in the first BZ as depicted by Fig. 0J All other parameters are the same as 
in Fig. |H1 




FIG. 10: Local "skeleton" diagrams for the DMFT self-energy S. Wavy lines represent the local 
(Hubbard) Coulomb interaction U, full lines denote the local Green function Gu. 




FIG. 11: Typical "skeleton" diagrams for the self-energy in the DMFT+Xk approach. The first 
two terms are DMFT self-energy diagrams; the middle two diagrams show contributions to the 
non-local part of the self-energy from spin fluctuations (see sectior JTTl|) represented as dashed lines; 
the last diagram (b) is an example of neglected diagram's leading to interference between the local 
and non-local parts. 
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FIG. 12: (Color online) Filling dependence of the pseudo gap potential A calculated with 
DMFT(QMC) for varying Coulomb interaction (U) and temperature (T) on a two-dimensional 
square lattice with two sets of (t,t'). 



